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We propose a scheme to mimic and directly measure the fractional particle number in a generalized 
Su-Schrieffer-Heeger model with ultracold fermions in one-dimensional optical lattices. We show that 
the fractional particle number in this model can be simulated in the momentum-time parameter 
space in terms of Berry curvature without a spatial domain wall. In this simulation, a hopping 
modulation is adiabatically tuned to form a kink-type configuration and the induced current plays 
the role of an analogous soliton distributing in the time domain, such that the mimicked fractional 
particle number is expressed by the particle transport. Two feasible experimental setups of optical 
lattices for realizing the required Su-Schrieffer-Heeger Hamiltonian with tunable parameters and 
time-varying hopping modulation are presented. We also show practical methods for measuring the 
particle transport in the proposed cold atom systems by numerically calculating the shift of the 
Wannier center and the center of mass of an atomic cloud. 
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I. INTRODUCTION 

Particle fractionalization has been recognized as a re¬ 
markable and fundamental phenomenon in both relativis¬ 
tic quantum field theory and condensed matter systems 
[Ml . The first physical demonstration of fractionaliza¬ 
tion is the celebrated Su-Schrieffer-Heeger (SSH) model 
of one-dimensional (ID) dimerized polymers @,0|, such 
as polyacetylene. In this model, a kink domain wall in 
the electron hopping configuration induces a zero-energy 
soliton state carrying a half-charge @ H]. The basic 
physics of fractionalization in SSH model is governed by 
a low-energy effective Dirac Hamiltonian with topologi¬ 
cally nontrivial background fields, which was firstly pro¬ 
posed by Jackiw and Rebbi [!, @. Subsequent achieve¬ 
ments were made to generalize the original SSH model 
to exhibit an irrational (arbitrary) fermion number by 
breaking the conjugation symmetry pjj-Hl. The frac¬ 
tional particle number (FPN) in these systems can be 
understood in terms of global deformations of the hole sea 
(or the valence band) due to the nontrivial background 
fields. 

SSH model has achieved great success in describing 
transport properties of polymers, and some novel phe¬ 
nomena associated with the topological solitons have also 
been explored in experiments 0. However, the intrin- 
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sic FPN has never been experimentally detected due to 
the spin-doubling problem in these materials: two spin 
orientations are present for each electron and thus a do¬ 
main wall in the polyacetylene carries an integer charge 
Q. Inspired by the newly discovered quantum spin Hall 
insulators 0, it was theoretically proposed to realize 
SSH model in an edge of this two-dimensional insulator 
by bringing a magnetic domain wall there, and the edge 
electrons with the inherent chiral symmetry may provide 
a direct signature of FPN 0 [l6j |. However, creation 
of such a magnetic domain wall acting only on the edge 
elections is experimentally challenging and the proposed 
schemes are yet to be demonstrated. 

In the past years, a lot of theoretical and experimental 
work has been carried out to simulate the Dirac equation 
and the involved exotic effects by using ultracold atoms 
PHI. Especially, it has been proposed to realize the 
(generalized) SSH model associated with effective Dirac 
Hamiltonian using ultracold atomic gases in the contin¬ 
uum 0 and in optical lattices [27H30j . The detection 
of FPN in these cold atom systems was also suggested 
by optical image of the density distribution of soliton 
modes 00. Since single-component fermionic gases 
or component-dependent optical lattices are used in the 
realization of FPN in atomic systems, the spin-doubling 
problem encountered in condensed matter systems can 
be avoided. In a recent experiment with a ID optical 
superlattice, SSH model in the absence of spatial domain 
walls was realized and its topological features were also 
probed [.31], making direct measurement of FPN in opti¬ 
cal lattices to be feasible and timely. 

In this work, we propose a new scheme to mimic and 
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directly measure FPN in the generalized SSH model us¬ 
ing ultracold fermions in ID optical lattices. Firstly, we 
show that FPN in this model can be simulated in the 
momentum-time parameter space in terms of Berry cur¬ 
vature without creating a spatial domain wall. In this 
simulation, a hopping modulation is adiabatically tuned 
to form a kink-type configuration and the induced cur¬ 
rent plays the role of an analogous soliton distributing in 
the time domain, so that the mimicked FPN in param¬ 
eter space is expressed by the adiabatic particle trans¬ 
port. Furthermore, we explore how to implement this 
new scheme with ultracold fermions in ID optical lat¬ 
tices. We propose two experimentally setups to realize 
the required SSH Hamiltonian with tunable parameters 
and hopping modulations, and then show practical meth¬ 
ods for measuring the particle transport in the proposed 
cold atom systems by numerically calculating the shift of 
the Wannier center and the center of mass of an atomic 
cloud. Some possible concerns in realistic experiments, 
such as the energy scales, the adiabatic condition and 
the effects of an external harmonic trap, are also consid¬ 
ered. In comparison with the previous proposals of re¬ 
alizing and detecting FPN [TH, l26l - l30j , one advantage of 
the presented scheme is that it does not require the spa¬ 
tial domain in the hopping configuration, which is usually 
hard to create. Another advantage is the adiabatic par¬ 
ticle transport corresponding to the value of PFN can be 
directly measured in our proposed cold atom systems. 

The rest of this paper is organized as follows: Section 
II presents a brief review on fractionlization in a gener¬ 
alized SSH model. In Sec. Ill, we elaborate our scheme 
of simulation and measurement of FPN in this model, 
based on the Berry curvature and adiabatic transport 
approaches. In section IV, we propose two feasible ex¬ 
perimental setups of ID optical lattices to realize the re¬ 
quired Hamiltonian, and then discuss how to measure 
the atomic particle transport in the proposed systems. 
Finally, a short conclusion is given in Sec. V. 

II. FRACTIONALIZATION IN SSH MODEL 

Before describing our scheme, we briefly review the ar¬ 
bitrary FPN in the generalized SSH model in this section. 
We start with this model described by a tight-binding 
Hamiltonian UM 

F = E[ J +(- 1 )^] (&Wi+H.c.) 

n 

+ ( 1 ) 

n 

where c n (cjj is the fermion annihilation (creation) op¬ 
erator in site n, J is the uniform hopping amplitude, <5 is 
the dimerized hopping modulation, and A is a staggered 
potential breaking the inversion symmetry [the conjuga¬ 
tion symmetry in the low-energy Dirac Hamiltonian (3)]. 
In this lattice system, as shown in Fig. 1(a), the even and 
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FIG. 1: (Color online) Two experimental setup for simulat¬ 
ing SSH Hamiltonian in ID optical lattices, (a) A double-well 
optical superlattice trapping noninteracting single-component 
fermionic atoms. A unit cell contains two nearest lattice 
sites with energy offset A (for atoms denoted by dark and 
light blue small balls) and the atomic hopping exhibits stag¬ 
gered modulation configuration, (b) A state-dependent opti¬ 
cal lattice trapping noninteracting two-component fermionic 
atoms, where atomic states | t) and | 1) are denoted by blue 
and red small balls. The block dotted line represents a Ra¬ 
man field Ri(x). (c) Raman-assisted tunneling. The uniform 
atomic hopping between the nearest-neighbour in the state- 
dependent optical lattice is realized by a large-detuned Raman 
transition with detuning Adi and Rabi frequencies of the Ra¬ 
man beams fii, 2 , while the hopping modulation is realized by 
another pair of Raman beams gu ,2 with large detuning Ad 2 - 
The Zeeman splitting gives rise to the A term in this system. 

odd number sites form two sublattices with modulating 
hopping amplitudes (on-site energies) and thus a unit cell 
contains two nearest lattice sites, which are used to con¬ 
stitute a pseudo-spin. When the inversion (conjugation) 
symmetry preserves with A = 0, the system corresponds 
to the original SSH model of polyacetylene. 

By employing Fourier transformation in the spin basis, 
we can obtain the Bloch Hamiltonian of the model as 

H = d(k) ■ <t, ( 2 ) 

where a = (a x ,a y ,a z ) are the Pauli matrices acting on 
the pseudo-spin, d(k) = (Jcos — <5sin A) is the 
band vector with a being the lattice spacing as shown in 
Fig. 1(a). By linearizing the Bloch bands near the Dirac 
point ko = tt/u, Hamiltonian (2) can be transformed into 
an effective low-energy relativistic Hamiltonian 

A 

Hd = VFPxCTx - 2S<Jy + —a z ( 3 ) 

in the continuum, where vf = Ja/h is the Fermi velocity, 
p x is the momentum operator measured from the Dirac 
point, <5 and A act as two background fields [Hi¬ 
lt has been widely studied that for a kink-type back¬ 
ground potential with S(x — > ±oo) = ±<5 0 , an unpaired 
soliton state appears at the kink carrying FPN [j, Til 

A/’s = — — arctan(-i^), (4) 

7r A 
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which may exhibit arbitrary fractional eigenvalues. The 
minus sign in FPN is due to the fact that the physical 
fermion number in the soliton sector is defined as be¬ 
ing measured relative to the free sector without the kink 
background, and this fractional part of the fermion num¬ 
ber actually comes from theglobal contribution (polar¬ 
ization) of the valence band [2. UlL l2a|. In addition, it has 
a topological character in the sense that it is dependent 
only on the asymptotic behavior of the background fields 
instead of their local profiles. When A —> 0, it recovers to 
the half fermion number ±1 for the zero-energy soliton 
mode in the original SSH model Q. 

It is interesting to note that fractionalization also ex¬ 
hibits in many low-dimensional correlated electron sys¬ 
tems. For instance, a well-known example is that of frac¬ 
tional excitations in the fractional quantum Hall regime 
1 , 0 , which is a consequence of strong Coulomb inter¬ 
action among 2D electrons in partially filled Landau lev¬ 
els. In addition, the collective excitations in some ID 
interacting fermion systems amay be characterized by 
effective fractional charges via the spin-charge separa¬ 
tion mechanism [10]. Fractionalization in these systems 
is basically due to election correlations and hence is com¬ 
pletely different from that in SSH model of noninteract¬ 
ing fermions. The fractional charges in some correlated 
electron systems have been directly observed in experi¬ 
ments [s!j, ■ However, despite the fractionalization in 

SSH model being investigated for decades, FPN there 
is yet to be directly measured in solid state materials 
(mostly due to the spin-doubling problem) or in artifi¬ 
cial systems, even for the simplest half-fernhon-number 
case with A = 0. Therefore, experimentally feasible 
schemes for direct measurement of the intrinsic FPN in 
SSH model would be of great value. 


III. SCHEME TO DIRECTLY MEASURE PFN 


Several schemes have been proposed to realize SSH 
model and to probe the soliton modes with FPN in cold 
atom systems [27H30| ; however, the realization of the re¬ 
quired kink background and the local detection of the 
fermion number of a soliton state therein are still chal¬ 
lenging in practical experiments. In the following, we 
propose a simpler scheme to mimic and then directly 
measure PFN via the adiabatic particle transport. 

The energy bands of SSH model described by Hamil¬ 
tonian (2) can be mapped onto a two-level system in 
the Bloch sphere with the parameterized Bloch vector 
S = ho (sin 9 cos <f>, sin0sin<(>, cosd), where 9 and <f> are 
respectively the polar and azimuthal angles. In this map¬ 
ping, we have 


ho = yj 4J 2 cos 2 (fca/2) + 4 5 2 sin 2 (fca/2) + A 2 /4, 

9 = arccos(A/2/i 0 ), (5) 

6 = — arctan 


S sin(ka/2 ) 
J cos(fca/2) 


The degeneracy point locates at k = n/a, 5 = 0 and A = 




FIG. 2: (Color online) (a) Berry curvature distribution of 
the ground band Tkt in the center region of the k-t parame¬ 
ter space. The Berry curvature outside is almost vanishing, 
and integration over the whole parameter space (0 ^ k ^ 
2n/a,0 ^ t ^ 2to) gives the particle transport, (b) Spin 
texture of SSH model, interpreted as a mapping from the k-t 
parameter space onto the Bloch sphere by Eq. (5). The colors 
show the SYcomponent cos[9(k, t)] and the arrows show the 
azimuthal component ( S x ,S y ) = (sin 9 cos (j), sin 6 sin 0). Here 
the lengths of all the arrows has been divided by a factor 7r 
for visibility, (c) The energy spectrum for a lattice system 
with size L = 100. (d) The variation of density in each lattice 
site with respect to the time-varying hopping modulation at 
hall filling, (e) The analog of kink potential and soliton-like 
current in the time domain under adiabatic conditions. The 
parameters in (a-e) are J = 1, <5o = A = 0.1 and to = 5/£. 


0. We consider the fractionalization in this band insula¬ 
tor at half filling [26|, corresponding to the low-energy 
level with the eigenstate |u_) = (sin |e - ^, — cos |) T 
with T being the transposition of matrix. In this 
framework, one can define the Berry connection Ag = 
(u-\idgu _) = 0 and A^ = (u-\idcf,u-) = sin 2 |. 

Instead of considering the kink background field in the 
spatial domain, here we introduce a time-varying hopping 
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modulation with a kink-type ramping configuration 

5(f) = 5 0 tanh[£(f-f 0 )], (6) 

where to denotes the center of the time domain wall and £ 
represents the ramp frequency. We assume to ^ l/£ (ac¬ 
tually to = 5/£ is large enough) such that 5(0) ~ —So at 
the beginning t = 0, and then adiabatically ramp the sys¬ 
tem to t = tf = 21 0 with S(tf) ~ 5 0 at the end. That is to 
say, we simulate a kink potential in the time domain in¬ 
stead of creating it in real space. In this dynamical case, 
the bulk gap is E g = 2-^/A 2 + 165(f) 2 , which will close 
at t = to for original SSH model with A = 0. To guaran¬ 
tee the adiabatic condition for original SSH case, which 
requires a gapped bulk band, we can use a time-varying 
staggered potential with the form A (t) = So sin(7rf/2f 0 ). 

We consider the adiabatic evolution of the system with 
the ramping parameter 5(f) (and A(f) for original SSH 
case) and provide that the Fermi level lies inside the band 
gap in the whole progression. The Berry phase effect of 
this ID band insulator can be measured from the par¬ 
ticle transport [j]. This is an analog of the adiabatic 
charge pumping proposed by Thouless [32j; however, the 
parametric driving in our case does not form a closed 
cycle but only a half one. The topological pumping in 
cold atom systems and photonic quasi-crvstals has been 
discussed in the contexts of SSH model |33l - [35l ] and ID 
quasi-periodic Harper model [36l - l40j , where the pumping 
particle is shown to be quantized one over one period and 
can be fractional over a fraction of one period [4Cj. 

In the momentum-time (k-t) parameter space, we can 
rewrite the Berry connection as A k = dkfyA^ + d k 0Ag 
and A t = dt^A^ + d t 0Ag. Thus the Berry curvature 
Fkt = dkA t — dtAk in the k-t space can be obtained as 


The particle transport for the ground band in this ID 
band insulator over the ramping progression of parameter 
5(f) from f = 0tof = f/is given by the integration of 
the Berry curvature 


Q = 




dkj~gt 


1 

27T 



dk [A k (k,tf) - A k (k, 0)] 


1 

=-arctan 

7r 
1 

cs-arctan 

7T 


4-7T J5o 


^A^tt 2 J 2 +45 2 + A 2 /4, 
45q' 


= K, 


( 8 ) 


where the approximation satisfies well for J 5o,A, 
and becomes exact when A = 0. In the calculation, the 
integration of d k A t over k vanishes due to the periodic 
condition in the Berry vector potential Q. Here the un¬ 
quantized adiabatic particle can be regarded as the po¬ 
larization change in this ID band insulator Hl|, which 
has also been discussed in the context of nanotubes and 
ferroelectrics materials [42j • If one consider an additional 
anti-kink-type modulation to form a full cycle, then the 
particle transport will be quantized (±1 or 0 depending 
on the loop of the cycle) after one period as the inte¬ 
gration of the Berry curvature in the extensional region 
contribute the other fractional portion jam. 

There are actually close connections between FPN of 
a soliton state and the adiabatic transport via the Berry 
phase approach in Eq. (8). We can consider the response 
equation in the progression of dynamical generation of 
the background field @,E1: 


T = _ AJsin 2 ^5(f) _ 

kt 2[4J 2 cos 2 !f + 45(f) 2 sin 2 % + A 2 /4]§ ’ 

where J and A are assumed to be constants here. We 
note that Berry curvature distribution given by Eq. (7) is 
modified for time-varying A (t) = So sm(nt/2to) discussed 
previously, however, the corresponding PFN given by fol¬ 
lowing equation (8) remains since it just depends on the 
boundaries of the background fields Q . 

Figure 2(a) shows an example of the Berry curvature 
distribution in the center region of the k-t space for typ¬ 
ical parameters, while the Berry curvature outside is 
almost vanishing. In small A limit, a sharp peak ex¬ 
hibits in the Berry curvature distribution at the posi¬ 
tion of ( k,t ) = (7r/a, to); which is the dominant contri¬ 
bution to its integration over the parameter space. Fig¬ 
ure 2(b) shows the corresponding spin texture, which is 
interpreted as a mapping from the k-t parameter space 
onto the Bloch sphere by Eq. (5). Figures 2(c) and 
2(d) show the static energy spectrum and the adiabatic 
density variation in each lattice site with respect to the 
hopping modulation for the lattice size L = 100, respec¬ 
tively. It looks like a density kink-soliton (anti-soliton) 
configuration appears in the time-domain. 


Ps = ^ d x Q{x,t ), j c = -^—d t Q(x, f), (9) 

where 0(cc,f) represents the generalized angular angle 
of the background field, p s denotes the soliton density 
distribution near the spatial domain wall, and j c is the 
induced current. In the present model with J> 5, A (in 
which case the same polarization variation is obtained by 
the band Hamiltonian and the low-energy Dirac Hamil¬ 
tonian), the angular angle 0 = — arctan(45/A) is just 
time-dependent. Therefore, the induced current mimics 
an analog of kink-soliton in the time domain, as shown 
in Fig. 2(e). The induced current over the whole time 
domain gives the transferred particle 

f tf 

Q= jc(t)dt , ( 10 ) 

J o 

which takes the value given by Eq. (8) and depends only 
on the boundaries of the angular angle under the adia¬ 
batic condition. 

So far, we have described our scheme to simulate FPN 
in SSH model in a parameter space and to directly mea¬ 
sure it via the adiabatic transport. In contrast to the 
previous schemes (TH], H(§j30|, the presented scheme does 
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not involve the spatial kink domain in the hopping con¬ 
figuration, which is usually hard to realize and (or) con¬ 
trol in experiments. In addition, the particle transport 
corresponding to the value of PFN can be directly mea¬ 
sured in cold atom systems, such as from the measure¬ 
ment of atomic density distribution and atomic current 
[43j ] . which will be discussed in ID optical lattices in the 
next section. 


IV. EXPERIMENTAL IMPLEMENTATION IN 
OPTICAL LATTICES 

In this section, we turn to discuss the implementation 
of our scheme of mimicking and measuring FPN in ID 
optical lattices. We first propose two experimental setups 
to realize the required SSH Hamiltonian with tunable pa¬ 
rameters, and then discuss how to measure the particle 
transport in the proposed cold atom systems by numeri¬ 
cally calculating the shift of the Wannier center and the 
center of mass of an atomic cloud. 


A. Two experimental setups 

The first experimental setup we proposed is a ID opti¬ 
cal supperlattice trapping a noninteracting atomic gas of 
single-component fermions, as shown in Fig. 1(a). Such 
an optical lattice has been widely realized in experiments 
[3ll 14314451 . It is generated by superimposing two lattice 
potentials with short and long wavelengths differing by a 
factor of two, with the optical potential given by 

V{x) = V\ sin 2 (fcix + y>) + V 2 sin 2 (2fci:r). (11) 


beams with a relative polarized angle [46| . The separa¬ 
tion and potential depth for different atomic components 
can be well controlled by the angle and the laser inten¬ 
sity, with a simple example of such a ID lattice potential 
as shown in Fig. 1(b) 

V a (x) = Vosin 2 (fc s a; ± 7 t/4). (12) 

Here Vo is the lattice potential depth, k s is the wavelength 
of the laser beams (the lattice spacing a = 2n/k s ), and 
±7 t/ 4 are the polarization angles for atomic states | t) 
and | |), respectively. 

For sufficiently deep lattices, the atoms in the system 
must alter their internal states in order to tunnel between 
two nearest-neighbor lattice sites. This can be achieved 
by the Raman-assisted tunneling method [H, |47| - |50| , as 
shown in Fig. 1(c). The energy offset between two atomic 
states arise from the external Zeeman field and then play 
the role of tunable parameter A in this system. Two 
pairs of Raman beams with Rabi frequencies Hi 2 and gi : 2 
are used to induce two large-detuned Raman transitions 
with detuning A^i and A d 2 , respectively. One can use 
the former pair of Raman beams to realize the uniform 
nearest-neighbor hopping 

J = A 0 J w^(x - x n )e lkxX wi(x - x n+1 )dx, (13) 

where Aq = IHiH^I/Adi is the effective Raman strength 
constant, k x is the momentum difference along the x axis 
between the two beams, and w a (x) are the Wannier func¬ 
tions of the lowest Bloch band for atomic state |cr). 

To realize the time-varying hopping modulation term, 
one can use another pair of laser beams with a resulting 
Raman field, as shown in Fig. 1(b), 


Here k\ is the wave vector of the short wavelength trap¬ 
ping lasers (the lattice spacing a = 2n/ki), ip and V\^ are 
respectively the relative phase and the strengths of the 
two standing waves. By varying the laser intensity and 
the phase, one can fully control the lattice system with 
ease [S3, l43H45j , and then make the system well described 
by Hamiltonian (1) of SSH model in the tight-binging 
regime (31) . In the experiments [U, I43i - l45j , the hopping 
configuration J + (— 1) J S can be adjusted by varying po¬ 
tential strengths Vj .,2 or swapping the relative phase p, 
and the staggered potential A can be tuned by the phase. 
Therefore in this system, a straightforward way to realize 
the required hopping modulation with kink-type config¬ 
uration in the time domain is by changing these tun¬ 
able parameters of the optical superlattice with a well- 
designed sequence 31]. 

Another experimental setup, which would be more con¬ 
venient as we will see in the following, is loading an ul¬ 
tracold Fermi gas of two-component (internal states |<r) 
with a =t, |) atoms in a state-dependent optical lattice 
[46| . It has been proposed to realize SSH Hamiltonian 
with a spatial domain wall in this system [27], and such 
a state-dependent optical lattice have been experimen¬ 
tally created by superposing two linearly polarized laser 


Ri =gigl/A d 2 = Ai(t)sin(2fc s a; + 7r/2), (14) 

where A\ is a time-dependent constant controlled by the 
laser intensities or detuning A d 2 - In this way, the hop¬ 
ping modulation (—l) n d(f) is given by 

A 1 (t) J w^(x — x n )sm(2k s x + n/2)w±(x — x n +i)dx. 

Here the staggered hopping modulation is a consequence 
of the relative spatial configuration of the lattice and the 
Raman field: the period of R\ ( x ) is double of the lattice 
period and R\{x) is antisymmetric corresponding to the 
center of each lattice site, as shown in Fig. 1(b). Thus, 
to realize the proposed particle transport scheme in this 
system, we can adjust the Zeeman field and the Raman 
field i? 0 to tune the parameter A and J in SSH Hamil¬ 
tonian, and then independently vary the intensity of an¬ 
other Raman field Ai(t) in time with a kink-type form. 
The case of time-varying A (t) can also be achieved in 
a similar way. The time modulation of the Raman cou¬ 
pling in ultracold atoms has been demonstrated in recent 
experiments [51] . 

Considering 40 K atoms and typical lattice spacing a = 
532 nm, one has the recoil energy Er/H ss 30 kHz. For 
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the optical superlattice system with intermediately deep 
lattice, a typical uniform hopping strength is J ~ 0.1 E R , 
and the other parameters <5 and A can be tunable in a 
wide regime [3l|. For the state-dependent optical lattice 
system, the uniform hopping strength given by Eq. (13) 
is proportional to the effective Raman intensity Ao, which 
is typically in order of megahertz, and the overlap integral 
of Wannier functions between neighbor lattice sites can 
be about 10 -2 [52j. Thus the Raman-induced uniform 
hopping strength in this system is J ~ 0AE R , and the 
nature (next nearest-neighbor) hopping fjv within sub¬ 
lattices can be effectively suppressed by sufficiently deep 
lattice Vo- For example, the numerical calculation shows 
that t n < 10 ~ 3 E r for V 0 « 22 E R [52|. We can consider 
typical parameters So and A (the minimum bulk gap in 
the dynamical progression is 2A) in the order of 0.1 J. In 
this case, the adiabatic approximation works well for the 
ramp time tf /i/0.1 J ~ 1 ms. Thus one can choose the 
ramp frequency £ = O.OlJ/h and to = 50 ms, which is 
well shorter than the typical coherence time in cold atom 
experiments. The non-adiabatic Landau-Zener transition 
from the ground band to the excited band for the chosen 
parameters is then given by Plz ~ / 4hS °£ ss 0. 

In addition, the finite temperature effects do not interfere 
with the particle transport progression for temperatures 
smaller than the energy gap [33j |. This requires the tem¬ 
perature of the order of 0.08 Er/Ub ~ 20 nk (ks is the 
Boltzmann constant), which has been achieved in current 
experiments with, e.g., 40 K atoms. So we can conclude 
that the required Hamiltonian with tunable parameters 
and the adiabatic condition are able to be realized under 
realistic circumstances. 


B. Experimental measurement methods 

It has been shown that the particle transport can be 
connected with the Wannier center based on the modern 
theory of charge polarization [3, E3 • Especially, the shift 
of the Wannier center in each unit cell is 


X d = x c (t = t f ) - x c (t = 0), (15) 

where x c = {w n \x\w n ) is the Wannier center, with 
| w n (x)) = A. e ~' lk ^ n ~ x '>\ii_(k)) being the Wannier 
function of the ground band in the n-th unit cell. The 
shift of the Wannier center encodes the adiabatic particle 
transport (the variation of polarization) as [13, HH, El 

X d /a = Q. (16) 

In Fig. 3, we have numerically calculated the variation 
of the Wannier center in each lattice site in the proposed 
system. For the symmetric case with A = 0 + (the solid 
blue line) and A (t) = So sin(7rf/2/ 0 ) (the dashed red line) 
shown in Fig. 3(a), after adiabatically tuning the hop¬ 
ping modulation over the kink form, the Wannier center 
in each lattice site shift downwards one site, that is, ex¬ 
actly one-half of the unit cell. According to Eqs. (8) 
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FIG. 3: (Color online) The shift of the Wannier center in each 
lattice as a response to the adiabatically tuning hopping mod¬ 
ulation. (a) The symmetric case with A = 0 + (solid blue line) 
and A(t) = do sin(7rt/2to) (dashed red line) with So = 0.1J. 
After tuning the hopping modulation over the kink form, the 
Wannier center for each lattice site shift downwards one site, 
i.e. one-half of the unit cell, corresponding to Q = —1/2 as 
expected for the half-charge in the original SSH model, (b) 
The symmetry-breaking case with A = 4<5o = 0.1J, the shift 
of the Wannier center is nearly half a lattice site, correspond¬ 
ing to Q = —1/4 expected in this case. Other parameters in 
(a) and (b) are J = 1 and to = 5/£. 


and (16), the adiabatic particle transport is Q = —1/2 
(the sign depends on the shift direction), as expected 
for the half-charge in the original SSH model. For the 
symmetry-breaking case with A = 4Sq = 0.1J as shown 
in Fig. 3(b), the shift of the Wannier center is nearly 
half a lattice site, which is consistent with the expected 
Q = —1/4 in this case. 

The shift of Wannier center shown in Fig. 3 implies the 
appearance of atomic current (the transport of particles) 
in each unit cell, which flow through the whole lattice sys¬ 
tem. Under the adiabatic condition, the atomic current 
will take the solitonic form in the ramping progression, 
which is similar to the example shown in Fig. 2(e). In 
principle, the transport dynamics can be detected by us¬ 
ing the sing le-atom in situ imaging technology in optical 
lattices [531 1 ■ Thus, the variation of atomic density dis¬ 
tribution in a unit cell associated with induced current 
can be experimentally extracted out in this way. In our 
proposed systems, it would be more convenient to detect 
the global current through the whole lattice by measur¬ 
ing the time evolution of the atom fractions in the even 
and odd sublattices, instead of using in situ detection in 
a single unit cell. For the double-well superlattice sys¬ 
tem in Fig. 1(a), the atomic current associated with 
the atom fractions of the even/odd sublattices has been 
measured in the experiment by transferring the atoms 
to higher-lying Bloch bands and applying a subsequent 
band mapping technique [43|. For the state-dependent 
optical lattice system in Fig. 1(b), the even and odd 
sublattices trap | t) and | 4-) atoms, respectively. There¬ 
fore in this system, one can simply measure the evolution 
of atom fractions via optical imaging the up-component 
(down-component) atoms, such as using state-resolved 
time-of-flight measurements [54|. Therefore, one can ob- 
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FIG. 4: (Color online) The shift of the center-of-mass of an 
atomic cloud as a function of (a) the parameter A with fixed 
Jo = 0.1; and (b) the parameter So with fixed A = 0.1, re¬ 
spectively, for a finite lattice system with L = 400 sites and at 
half filling N a /L = 1/2. The solid (blue) line and the dashed 
(black) line denotes the corresponding particle transport Q 
obtained by Eq. (8) and fractional particle number A f s ob¬ 
tained by Eq. (4) in the text, respectively. Other parameters 
in (a) and (b) are J = 1 and to = 5/£. 


tain the particle transport Q by its integration over the 
time domain, as given by Eq. (10), which corresponding 
to the mimicked FPN in this system. 

In addition, the particle transport can be directly mea¬ 
sured from the shift of the center of mass of an atomic 
cloud in a finite lattice system (diil. [35l[40| . In our system, 
the center of mass of an atomic cloud in the lattice with 
L sites x com (t) is given by 

1 L 

%com (t) = ^ ^ ^ ^ t ) | 71, (17) 

... — i - 


where N a = L/ 2 is the atomic number at half filling, e oc 
denotes the occupied state of the fermionic atoms, and 
tpoc is the corresponding wave function. Under the adia¬ 
batic evolution with J(0) —>• S(tf) described by Eq. (6), 
the center of mass of the system shift from the position 
Xcom(0) to i com (tf). It can be proved that the shift of 
the center of mass 5x c om = x com (tf) — a.’ CO m(0) is pro¬ 
portional to the particle transport in the infinite L limit 

mm 


yfacom = Q■ (18) 

If L is large enough, such that the bulk properties of the 
system are almost not affected by the edges, N a 5x coln /L 
in the above equation will be approximated to be the 
ideal particle transport Q in infinite system. In Fig. 4, 
we have calculated the shift of the center of mass of an 
atomic cloud for a lattice system with L = 400 sites, 
N a Sx C om/L (the red circles) as a function of the param¬ 
eters A [Fig.4 (a)] and Jo [Fig-4 (b)], respectively. As 
shown in Fig. 4(a) and 4(b), the calculated shift of the 
center of mass is well described by the particle transport 
Q (the solid blue line) obtained by Eq. (8) with small 
deviations. These deviations are due to the finite size 


effects and become smaller and smaller with the increase 
of the lattice size in our simulations, which will be invis¬ 
ible in practical experiments. For comparison, we also 
plot the mimicked FPN J\f s given by Eq. (4) in this sys¬ 
tem in Fig. 4 (the dashed black line). From Fig. 4, 
one can see that the corresponding FPN is nearly equal 
to the particle transport within the parameter regimes. 
In current experiments, the center-of-mass position of an 
atomic cloud can be directly and precisely measured, ei¬ 
ther by using in situ measurement of the atomic density 
distribution in the optical lattices [HJ or deduced from 
the time-of-flight imaging f54j . 

Finally, we note that a shallow-enough harmonic trap 
in practical experiments will not affect the particle trans¬ 
port [33|,[39j and the main results of this paper remain in¬ 
tact. In order to take the effect of the harmonic trap into 
account, we can add a term H t = V t ~ L/2) 2 cl l c n 

into Hamiltonian (1), where V) is the trap strength and L 
is the lattice size. Within a local density approximation, 
the lower band will still be filled at the center of the trap 
and thus the shift of the Wannier center in these lattice 
sites is expected to be nearly the same as those shown 
in Fig. 3; while the band is only partially filled near the 
edge with the local trap energy V t (n — L/2) 2 > E g , such 
that the pumping argument does not apply to this region 
[33-j. Therefore in practical experiments, one may empha¬ 
size on the shift of Wannier center in the central region 
or turn the trap strength to a small value V) ~ AE g /L 2 . 
For the shift of the center of mass of an atomic cloud, 
our numerical simulations demonstrate that the results 
shown in Fig. 4 preserve with a deviation less than 2% 
for V t ss 0.6 x 10~ 5 J, while V t ~ 10~ 4 J is enough if the 
lattice size reduces to L = 100, which are consistent with 
the estimates in the local-density analysis. 


V. CONCLUSIONS 

In summary, we have proposed a scheme to mimic and 
measure the FPN in the generalized SSH model with cold 
fermions in ID optical lattices. It has been shown that 
FPN in this model can be simulated in the momentum¬ 
time parameter space in terms of Berry curvature with¬ 
out a spatial domain wall. In this simulation, a hopping 
modulation is adiabatically tuned to form a kink-type 
configuration and the induced current plays the role of 
soliton in the time domain, so FPN is expressed by the 
particle transport. We have also proposed two experi¬ 
mental setups of optical lattices to realize the required 
Hamiltonian with tunable and time-varying parameters, 
and considered the energy scales and the adiabatic con¬ 
dition under practical circumstances. Finally we have 
discussed how to directly measure the particle transport 
in the proposed systems by numerically calculating the 
shift of the Wannier center and the center of mass of an 
atomic cloud. Considering that all the ingredients to im¬ 
plement our scheme in optical lattices have been achieved 
in the recent experiments, it is anticipated that the pre- 















sented proposal will be tested in an experiment in the 
near future. The direct measurement of such a mimicked 
FPN in cold atom experiments will be an important step 
toward exploring fractionalization and topological states 
in cold atom systems. Extensions of this work can en¬ 
able to simulate and measure FPN emerging in two- and 
three-dimensional Dirac Hamiltonian with topologically 
nontrivial (vortex and monopole) background fields [lj, [8j, 
which has been theoretically studied but remains elusive 
in nature. It will be also interesting to simulate a variety 
of topological states @] and study their properties in the 


parameter space using cold atoms. 


VI. ACKNOWLEDGEMENTS 


This work was supported by the RGC of Hong Kong 
(Grants No. HKU7045/13P and HKU173051/14P), the 
SKPBR of China (Grants No. 2011CB922104 and 
2013CB921804), the NSFC (Grants No. 11125417 and 
11474153), and the PCSIRT (Grant No. IRT1243). 


[1] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 

[2] A. Niemi and G. Semenoff, Phys. Rep. 135, 99 (1986). 

[3] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. 
Lett. 42, 1698 (1979); W. P. Su, J. R. Schrieffer, and A. 

J. Heeger, Phys. Rev. B 22, 2099 (1980). 

[4] A. J. Heeger, S. Kivelson, J. R. Schrieffer, and W. P. Su, 
Rev. Mod. Phys. 60, 781 (1988). 

[5] R. B. Laughin, H. Stormer, and D. Tsui, Rev. Mod. Phys. 
71, 863 (1988). 

[6] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
78, 195424 (2008). 

[7] D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 
1959 (2010). 

[8] C.-Y. Hou, C. Chamon, and C. Mudry, Phys. Rev. Lett. 
98, 186809 (2007); B. Seradjeh, C. Weeks, and M. Franz, 
Phys. Rev. B 77, 033104 (2008); A. Riiegg and G. A. 
Fiete, Phys. Rev. B 83, 165118 (2011). 

[9] R. de-Picciotto, M. Reznikov, M. Heiblum, V. Umansky, 
G. Bunin, and D. Mahalu, Nature (London) 389, 162 
(1997). 

[10] H. Steinberg, G. Barak, A. Yacoby, L. N. Pfeiffer, K. W. 
West, B. I. Halperin, and K. L. Hur, Nat. Phys. 4, 116 
(2008); H. Kamata, N. Kumada, M. Hashisaka, K. Mu- 
raki, and T. Fujisawa, Nat. Nanotechnol. 9, 177 (2014). 

[11] J. Goldstone and F. Wilczek, Pliys. Rev. Lett. 47, 986 
(1981); H. Takayama, Y. R. Lin-Liu, and K. Maki, Phys. 
Rev. B 21, 2388 (1980). 

[12] M. J. Rice and E. J. Mele, Phys. Rev. Lett. 49, 1455 
(1982). 

[13] R. Jackiw and G. Semenoff, Phys. Rev. Lett. 50, 439 
(1983). 

[14] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010); X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 
1057 (2011). 

[15] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Nat. Phys. 4, 
273 (2008). 

[16] J. I. Vayrynen and T. Ojanen, Phys. Rev. Lett. 107, 
166804 (2011). 

[17] S.-L. Zhu, B. Wang, and L.-M. Duan, Phys. Rev. Lett. 
98, 260402 (2007); C. Wu, Phys. Rev. Lett. 101, 186807 
(2008); L.-K. Lim, C. M. Smith, and A. Hemmerich, 
Phys. Rev. Lett. 100, 130402 (2008); J.-M. Hou, W.-X. 
Yang, and X.-J. Liu, Phys. Rev. A 79, 043621 (2009); 

K. L. Lee, B. Gremaud, R. Han, B. G. Englert, and 
C. Miniatura, Phys. Rev. A 80, 043411 (2009); A. 
Bermudez, N. Goldman, A. Kubasiak, M. Lewenstein, 
M. A. Martin-Delgado, New J. Phys. 12, 033041 (2010). 

[18] N. Goldman, A. Kubasiak, A. Bermudez, P. Gaspard, M. 


Lewenstein, and M. A. Martin-Delgado, Phys. Rev. Lett. 
103, 035301 (2009); A. Bermudez, L. Mazza, M. Rizzi, 
N. Goldman, M. Lewenstein, and M. A. Martin-Delgado, 
Phys. Rev. Lett. 105, 190404 (2010). 

[19] J. Y. Vaishnav and C. W. Clark, Phys. Rev. Lett. 100, 
153002 (2008); G. Juzeliunas, J. Ruseckas, M. Lindberg, 

L. Santos, and P. Ohberg , Phys. Rev. A 77, 011802(R) 
(2008); Y. Zhang, L. Mao, and C. Zhang, Phys. Rev. 
Lett. 108, 035302 (2012); S. Cao, D.-W. Zhang, H. Yan, 
and Z.-Y. Xue, JETP Lett. 98, 858 (2014). 

[20] M. Merkl, A. Jacob, F. E. Zimmer, P. Ohberg, and L. 
Santos, Phys. Rev. Lett. 104, 073603 (2010); S.-L. Zhu, 
D.-W. Zhang, and Z. D. Wang, Phys. Rev. Lett. 102, 
210403 (2009); D.-W. Zhang, Z.-Y. Xue, H. Yan, Z. D. 
Wang, and S.-L. Zhu, Phys. Rev. A 85, 013628 (2012). 

[21] See the reviews, D.-W. Zhang, Z. D. Wang, and S.-L. 
Zhu, Front. Phys. 7, 31 (2012); F. Mei, D.-W. Zhang, 
and S.-L. Zhu, Chin. Phys. B 22, 116106 (2013). 

[22] T. Salger, C. Grossert, S. Kling, and M. Weitz, Phys. 
Rev. Lett. 107, 240401 (2011). 

[23] L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu, T. 
Esslinger, Nature (Landon), 483, 303 (2012). 

[24] L. J. LcBlanc, M. C. Beeler, K. Jimenez-Garcia, A. R. 
Perry, S. Sugawa, R. A. Williams, and I. B. Spielman, 
New J. Phys. 15, 073011 (2013). 

[25] C. Qu, C. Hamner, M. Gong, C. Zhang, and P. Engels, 
Phys. Rev. A 88, 021604 (2013). 

[26] D.-W Zhang, L.-B Shao, Z.-Y Xue, H. Yan, Z. D. Wang, 
and S.-L. Zhu, Phys. Rev. A 86, 063616 (2012). 

[27] J. Ruostekoski, G. V. Dunne, and J. Javanainen, Phys. 
Rev. Lett. 88, 180401 (2002); J. Ruostekoski, J. Ja¬ 
vanainen, and G. V. Dunne, Phys. Rev. A 77, 013603 
(2008). 

[28] X.-J. Liu, Z.-X. Liu, and M. Cheng, Phys. Rev. Lett. 
110, 076401 (2013). 

[29] X. Li, E. Zhao, and W. V. Liu, Nat. Commun. 4, 1523 
(2013). 

[30] J. Javanainen and J. Ruostekoski, Phys. Rev. Lett. 91, 
150404 (2003). 

[31] M. Atala, M. Aidelsburger, J. T. Barreiro, D. Abanin, 
T. Kitagawa, E. Dernier, and I. Bloch, Nat. Phys. 9, 795 
(2013). 

[32] D. J. Thouless, Phys. Rev. B 27, 6083 (1983). 

[33] L. Wang, M. Troyer, and X. Dai, Phys. Rev. Lett. Ill, 
026802 (2013). 

[34] H. Guo, Phys. Lett. A 378, 1316 (2014); H. Guo, arXiv: 
1311.4134. 



9 


[35] Y. Qian, M. Gong, and C. Zhang, Phys. Rev. A 84, 
013608 (2011). 

[36] P. G. Harper, Proceedings of the Physical Society. Sec¬ 
tion A 68, 874 (1955); S. Aubry and G. Andre, Ann. Isr. 
Phys. Soc. 3, 133 (1980). 

[37] Y. E. Kraus, Y. Lahini, Z. Ringel, M. Verbin, and O. 
Zilberberg, Phys. Rev. Lett. 109, 106402 (2012); Y. E. 
Kraus and O. Zilberberg, Phys. Rev. Lett. 109, 116404 
(2012); M. Verbin, O. Zilberberg, Y. E. Kraus, Y. Lahini, 
and Y. Silberberg, Phys. Rev. Lett. 110, 076403 (2013). 

[38] L.-.J. Lang, X. Cai, and S. Chen, Phys. Rev. Lett. 108, 
220401 (2012); F. Mei, S.-L. Zhu, Z. M. Zhang, C. H. Oh, 
and N. Goldman, Phys. Rev. A 85, 013638 (2012). 

[39] F. Mei, J.-B. You, D.-W. Zhang, X. C. Yang, R. Fazio, 
S.-L. Zhu, and L. C. Kwek, Phys. Rev. A 90, 063638 
(2014). 

[40] P. Marra, R. Citro, and C. Ortix, Phys. Rev. B 91, 
125411 (2015). 

[41] R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 
1651 (1993). 

[42] E. J. Mele and P. Krai, Phys. Rev. Lett. 88, 056803 
(2002); S. Onoda, S. Murakami, and N. Nagaosa, Phys. 
Rev. Lett. 93, 167602 (2004). 

[43] M. Atala, M. Aidelsburger, M. Lohse, J. T. Barreiro, 
B. Paredes, and I. Bloch, Nat. Phys. 10, 588 (2014); 
M. Aidelsburger, M. Atala, M. Lohse, J. T. Barreiro, 
B. Paredes, and I. Bloch, Phys. Rev. Lett. Ill, 185301 
(2013); M. Aidelsburger, M. Atala, S. Nascimbene, S. 
Trotzky, Y.-A. Chen, and I. Bloch, Phys. Rev. Lett. 107, 
255301 (2011). 

[44] J. Sebby-Strabley, M. Anderlini, P. S. Jessen, and J. 
V. Porto, Phys. Rev. A 73, 033605 (2006); J. Sebby- 
Strabley, B. L. Brown, M. Anderlini, P. J. Lee, W. D. 
Phillips, J. V. Porto, and P. R. Johnson, Phys. Rev. Lett. 
98, 200405 (2007). 

[45] S. Foiling, S. Trotzky, P. Cheinet, M. Feld, R. Saers, A. 
Widera, T.Muller, and I. Bloch, Nature (London) 448, 
1029 (2007); S. Trotzky, P. Cheinet, S. Foiling, M. Feld, 
U. Schnorrberger, A. M. Rey, A. Polkovnikov, E. A. Dern¬ 
ier, M. D. Lukin, and I. Bloch, Science 319, 295 (2008). 


[46] O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. Han- 
sch, and I. Bloch, Phys. Rev. Lett. 91, 010407 (2003); P. 
J. Lee, M. Anderlini, B. L. Brown, J. Sebby-Strabley, W. 
D. Phillips, and J. V. Porto, Phys. Rev. Lett. 99, 020402 
(2007). 

[47] D. Jaksch and P. Zollcr, New J. Phys. 5, 56 (2003); E. J. 
Mueller, Phys. Rev. A 70, 041603 (2004). 

[48] F. Gerbier and J. Dalibard, New J. Phys. 12, 033007 
(2010); L. Mazza, A. Bermudez, N. Goldman, M. R.izzi, 
M. Angel Martin-Delgado, and M. Lewenstein, New J. 
Phys. 14 015007 (2012). 

[49] H. Miyake, G. A. Siviloglou, C. J. Kennedy, W. C. Bur¬ 
ton, and W. Ketterle, Phys. Rev. Lett. Ill, 185302 
(2013). 

[50] See the reviews, J. Dalibard, F. Gerbier, G. Juzeliunas, 
and P. Ohberg, Rev. Mod. Phys. 83, 1523 (2011); V. 
Galitski and I. B. Spielman, Nature (London) 494, 49 
(2013); N. Goldman, G. Juzeliunas, P. Ohberg, and I. B. 
Spielman, Rep. Prog. Phys. 77, 126401 (2014); H. Zhai, 
Rep. Prog. Phys. 78, 026001 (2015); S. Zhang, W. S. 
Cole, A. Paramekanti, and N. Trivedi, arXiv: 1411.2297. 

[51] K. Jimenez-Garcia, L. J. LeBlanc, R. A. Williams, M. 
C. Beeler, C. Qu, M. Gong, C. Zhang, and I. B. Spiel¬ 
man, Phys. Rev. Lett. 114, 125301 (2015); Y.-J. Lin, R. 
L. Compton, K. Jimenez-Garcia, W. D. Phillips, J. V. 
Porto, and I. B. Spielman, Nat. Phys. 7, 531 (2011). 

[52] J. K. Pachos, E. Alba, V. Lahtinen, and J. J. Garcia- 
Ripoll, Phys. Rev. A 88, 013622 (2013); D.-W. Zhang, 
C.-J. Shan, F. Mei, M. Yang, R.-Q. Wang, and S.-L. Zhu, 
Phys. Rev. A 89, 015601 (2014). 

[53] N. Gemelke, X. Zhang, C.-L. Hung, and C. Chin, Nature 
(London) 460, 995 (2009); W. S. Bakr, J. I. Gillen, A. 
Peng, S. Foiling, and M. Greiner, Nature (London) 462, 
74 (2009). 

[54] C. D. Fertig, K. M. OHara, J. H. Huckans, S. L. Rolston, 
W. D. Phillips, and J. V. Porto, Phys. Rev. Lett. 94, 
120403 (2005); N. Strohmaier, Y. Takasu, K. Gunter, R. 
.Jordens, M. Kohl, H. Moritz, and T. Esslinger, Phys. 
Rev. Lett. 99, 220601 (2007). 



